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ABSTRACT 

Driven  by  the  need  for  identifying  the  presence  of  flow  noise 
components  in  a  towed  volumetric  acoustic  array,  methodolo¬ 
gies  are  developed  for  high  resolution  3 -dimensional  (3-D) 
wavenumber-frequency  (k-f)  analysis  based  on  Krylov  space 
techniques.  These  methods  are  applied  to  actual  noise  data 
from  a  sea  trial.  The  experimental  results  indicate  that  the 
Krylov-based  method  outperforms  conventional  approaches 
and  provides  efficient  and  effective  solutions  to  k-f  analysis. 

Index  Terms —  Sonar  Signal  Processing,  Array  Signal 
Processing,  Sonar  Imaging,  Krylov  Space. 

1.  INTRODUCTION  AND  BACKGROUND 

In  order  to  reduce  the  interfering  effects  of  noise  radiated  by 
the  platform  ship,  passive  and  active  sonar  systems  gener¬ 
ally  utilize  arrays  of  hydrophones  mounted  on  a  cable  that 
are  towed  behind  the  ship  as  their  receivers.  Besides  water¬ 
borne  acoustic  noise  (e.g.  ambient  noise,  tow  ship,  shipping, 
biologies  etc.),  another  important  source  of  noise  or  interfer¬ 
ence  present  in  towed  arrays  is  flow  noise.  Flow  noise  is  non¬ 
acoustic  in  nature  and  is  caused  by  boundary  layer  turbulence 
exciting  the  array  structure  (including  excitation  of  flexural- 
like  waves  in  the  array  hose  wall)  and  inducing  pressure  fluc¬ 
tuations  on  the  surface  of  the  hydrophones  [1,2].  Flow  noise 
can  significantly  degrade  sonar  performance,  e.g.  probabil¬ 
ity  of  target  signal  detection,  especially  at  high  tow  speeds. 
Therefore  its  presence  and  properties  are  of  much  interest  to 
sonar  system  designers. 

This  work  is  motivated  by  the  problem  of  identifying  the 
presence  of  flow  noise  in  a  small  sparsely-sampled  volumet¬ 
ric  array.  Much  insight  into  the  types  of  noise  present  in  an 
acoustic  array  can  be  obtained  by  performing  a  wavenumber- 
frequency  (i k-f)  analysis  of  the  space-time  noise  field  mea¬ 
sured  by  the  array  as  it  is  towed  through  the  water.  A  volu¬ 
metric  array  requires  3-D  wavenumber  analysis  as  a  function 
of  frequency  because  of  the  3-D  arrangement  of  hydrophones. 
Hence  we  begin  by  reviewing  3-D  k-f  analysis. 

Let  P(pm,  t)  denote  the  signal  (acoustic  and  flow  noise) 
time  series  measured  at  time  t  by  the  ra-th  hydrophone  lo¬ 
cated  at  position  pm.  Then  the  wavenumber- frequency  power 
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spectrum  is  defined  as  the  expected  value 


p(k,/)  =  sj 

Or  equivalently, 


Y  [  P(Pm,t)e~^kTP--^dt 

m  ’’ 

P(k,/)  =  e"(k)R[/]e(k), 


(1) 

(2) 


where  eH(k)  =  [e^^P1  eJ27rkTp2  •  •  •  e'27rkTp"  ]  is  a  3-D 

wavenumber  steering  vector;  vector  p m  contains  3-D  posi¬ 
tion  coordinates  of  the  m-th  array  element,  and  R[f]  is  the 
nominal  spatial  covariance  matrix  of  the  noise  field  measured 
by  the  array  at  frequency  /.  In  the  above  equations,  the  3-D 
wavenumber  vector  k  =  [kx  ky  kz]T  is  defined  in  a  Cartesian 
coordinate  system.  In  this  paper,  kx  is  always  measured  down 
the  array  axis.  From  formula  (2)  we  see  that  the  wavenumber- 
frequency  spectrum  is  nothing  more  than  a  three-dimensional 
spatial  Fourier  analysis  of  the  measured  hydrophone  signa¬ 
tures  as  a  function  of  temporal  frequency. 

The  structure  of  the  wavenumber-frequency  power  spec¬ 
trum  is  highly  revealing  about  the  types  of  noise  sources  present. 
This  is  best  illustrated  using  a  equi-spaced  line  array  example 
which  needs  only  a  1-D  k-f  analysis,  that  is,  as  a  function  of 
kx  and  /.  When  the  array  elements  are  separated  in  distance 
by  d,  all  plane  wave  acoustic  signals  must  all  lie  in  the  k-f 
cone  depicted  in  fig.  1.  The  region  within  the  cone  is  known 
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Fig.  1.  Idealized  k  —  f  spectrum:  acoustic  and  non-acoustic  regions. 

as  the  visible  acoustic  region.  For  example,  an  acoustic  plane 
wave  signal  impinging  upon  the  array  from  direction  6  will 
manifest  itself  as  a  ridge-like  structure  that  lies  within  the 
acoustic  cone  as  shown  in  fig.  1.  Far  field  acoustic  noise  also 
lies  within  the  acoustic  cone  or  visible  region.  However,  non¬ 
acoustic  noise  components  such  as  flow  noise  will  be  present 
outside  the  acoustic  cone  and  often  has  a  ridge-like  structure. 
Hence  the  presence  of  non-acoustic  noise  in  an  array  can  often 
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be  quickly  determined  by  performing  a  simple  wavenumber- 
frequency  analysis.  Similar  results  apply  to  a  volumetric  ar¬ 
ray  except  that  the  visible  acoustic  region  now  forms  a  4-D 
hypercone  in  the  space  defined  by  the  coordinates  kx,  ky,  kz, 
and/. 

Estimating  the  k-f  spectrum  from  real  at-sea  data  is  hard. 
Real  acoustic  data  is  often  contaminated  by  acoustic  inter¬ 
ference  such  as  that  arising  from  the  tow  ship  (as  a  result, 
sidelobe  leakage  from  interferers  can  obscure  weaker  k-f  fea¬ 
tures)  and  can  be  non- stationary  because  of  array  motion.  In 
the  case  of  a  volumetric  array,  an  especially  challenging  prob¬ 
lem  is  visualization  of  the  4-D  k-f  volume  on  a  2-D  display. 

In  the  remainder  of  the  paper,  we  develop  a  Krylov  space- 
based  high  resolution  k-f  spectrum  estimator  and  apply  it  to 
actual  at-sea  data  collected  from  a  volumetric  array  to  identify 
the  presence  of  flow  noise.  The  Krylov-based  high  resolution 
estimator  is  compared  against  the  periodogram  and  Capon’s 
method  [4]  and  is  shown  to  work  well.  We  also  discuss  and 
demonstrate  methods  for  visualizing  the  4-D  k-f  volume  using 
slices  of  the  spectrum  taken  in  a  spherical  coordinate  system. 

2.  PROPOSED  SOLUTIONS  AND  RESULTS 

The  array  of  interest  is  comprised  of  triplet  sets  of  hydrophones 
that  are  arranged  and  mounted  in  a  triangular-like  pattern  (see 
fig.  2)  on  semi-rigid  supports  that  are  embedded  in  a  hose¬ 
like  structure.  Unlike  a  conventional  line  array,  the  volumetric 
arrangement  of  hydrophones  permits  unambiguous  left-right 
identification  of  acoustic  sources. 


y 


(b)  - 1 - ; - y 

Fig.  2.  Schematic  diagram  of  the  volumetric  array;  (a)  Full  view, 
(b)  View  down  the  array  axis. 

The  conventional  periodogram-based  method  generates  a 
4-D  (3-D  in  space  or  wavenumber  and  1-D  in  temporal  fre¬ 
quency)  volume  revealing  the  wavenumber-frequency  distri¬ 
bution  of  signal  modes.  That  is, 

PPER(k,f,t)  =  eH(k)R[f,t]e(k).  (3) 

Different  from  eq.(2),  here  R[/,  t\  is  a  time-varying  sample 
spatial  covariance  (the  corss-spectral  density  matrix  CSDM) 
obtained  from  a  finite  number  of  array  noise  measurements. 
It  will  be  helpful  for  the  purposes  of  visualizing  the  4-D  k-f 
volume  to  represent  the  wavenumber  vector  defined  in  eq.  (1) 
using  spherical  coordinates  as  k  =  k  u(</>,  0)  where  u(</>,  6)  = 
[sin  9  cos  4>  sin  6  sin  </>  cos  0] T  is  the  directional  cosine  vec¬ 
tor  as  a  function  of  the  angles  6  and  </>.  When  the  physical 
wavenumber  ||k||F  =  k  lies  within  the  visible  acoustic  re¬ 
gion  at  a  given  frequency,  the  angles  cf)  and  0  correspond  to 


the  actual  azimuthal  and  elevation  angles  respectively  of  an 
acoustic  signal  impinging  upon  the  array. 

High  levels  of  side-lobe  leakage  and  limited  resolution 
power  associated  with  periodogram  approach  affects  the  per¬ 
formance  of  wavenumber- frequency  analysis  in  revealing  use¬ 
ful  spatial  modes  and  flow  noise,  especially  when  strong  tow- 
ship  interference  is  present  in  the  estimated  CSDM  R [/,£]. 
Fig.  3(a)  shows  a  2-D  slice  cut  out  of  the  4-D  k-f  volume  in  (3) 
obtained  by  setting  (0,  <f)  =  (90°,  0°)  and  plotting  wavenum¬ 
ber  k  versus  frequency  /  in  the  spherical  coordinate  system 
(this  is  equivalent  to  scanning  kx  versus  frequency  while  set¬ 
ting  ky  —  kz  =  0).  The  strong  ridge-like  feature  present  in 
the  k-f  spectrum  that  roughly  follows  the  lower  dashed  line  is 
radiated  noise  from  the  tow  ship.  The  weaker  ridges  in  the 
same  orientation  in  the  acoustic  and  non-acoustic  regions  are 
a  result  of  sidelobes.  This  clearly  illustrates  the  deleterious 
effects  of  side-lobe  leakage  of  periodogram-like  approaches. 

Greatly  improved  performance  in  terms  of  sidelobe  leak¬ 
age  and  resolution  can  be  obtained  by  using  a  high  resolu¬ 
tion  spectral  estimation  method  such  as  Capon’s  method  [4]. 
Capon’s  method  can  be  directly  applied  to  estimate  the  4-D 
wavenumber- frequency  volume, 


,(k  ,/,*) 


1 

eH(k)  R-1  [/,  t\  e(k)  * 


(4) 


However,  inverting  the  R [/,  t\  with  a  large  eigen-spread  will 
result  in  mode  smearing  and  focusing  problems.  Shown  in 
Fig.  3(b)  is  a  2-D  slice  of  k-f  spectrum  from  the  4-D  k-f  analy¬ 
sis  results  ((9,  </>)  =  (90°,  0°))  using  Capon’s  method  in  eq. 
(4).  Although  greatly  improved  sidelobe  leakage  over  the  pe¬ 
riodogram,  we  will  show  in  the  next  section  that  even  bet¬ 
ter  performance  can  be  attained  using  vector  conjugate  gradi¬ 
ent  (V-CG)  rank-reduction  approaches. 


2.1.  Improving  performance  using  V-CG  rank-reduction 


In  [8],  motivated  by  the  problem  of  adaptively  detecting  a  sig¬ 
nal  using  an  array  with  a  large  number  of  sensors,  we  advo¬ 
cated  using  a  Krylov  subspace-based  adaptive  dimensionality 
reduction  technique  to  improve  detection  performance  and  re¬ 
duce  sensitivity  to  model  order  uncertainty.  In  this  applica¬ 
tion  we  are  using  the  array  manifold  vector  to  scan  the  3-D 
wavenumber  space  to  reveal  interesting  activities.  Roughly, 
only  those  noise  and  interference  spatial  modes  (both  acoustic 
and  non-acoustic)  that  are  not  orthogonal  to  the  steering  replica 
or  array  manifold  are  relevant  to  estimating  the  spectrum  at 
that  direction.  Such  observation  motivates  us  to  propose  and 
study  an  effective  procedure  that  directly  formulates  a  sub¬ 
space  for  capturing  only  the  spatial  modes  relevant  to  a  given 
steering  vector  pointing  to  a  direction  of  interest,  i.e.,  rank 
reduction.  Therefore,  to  alleviate  resolution  problems  caused 
by  inverting  a  high-dimensional  CSDM  R [/,  t\  with  a  large 
eigen-spread,  we  propose  using  rank  reduction  based  on  the 
V-CG  method.  One  then  projects  the  estimated  CSDM  R [/,  t\ 
onto  the  V-CG  estimated  subspace  to  yield  an  improved  struc¬ 
tural  estimate  of  the  CSDM  for  wavenumber- frequency  analy- 
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sis.  Therefore,  we  can  directly  use  the  V-CG  method  in  for¬ 
mulating  the  projected  structural  estimate  of  CSDM,  yet  at 
the  same  time  to  control  the  rank  of  the  space-time  V-CG  fil¬ 
ter  to  avoid  troubles  caused  by  the  large  eigen-spreads  of  the 
sample  estimate  R [/,£].  Hence,  a  sequence  of  4-D  rank-r 
wavenumber- frequency  spectra  can  be  obtained  using  a  rank- 
r  V-CG  method, 


}( rank—r ) 


(•<•/./)  = 


1 


(5) 


erf(k)  wr(k,  f,  t)’ 
where  the  rank-r  V-CG  space-time  filter  vector  wr(k,  /,  t)  is 
constructed  out  of  the  Krylov  space  according  to 


r 


wr(k, /, t)  =  G  £r{R[/,t],e(k)},  (6) 

i=l 

with  Krylov  space  JCr{ R,  e}  =<  e,  Re,  •  •  •  ,  Rr  1  e  >.  In 
the  above  V-CG  space-time  filter  formulation,  the  conjugate 
direction  vectors  d*  and  combination  coefficients  ai  are  adap¬ 
tively  calculated  from  the  estimated  CSDM  and  the  array  man¬ 
ifold  specified  steering  vector,  without  matrix  inversion,  using 
an  iterative  procedure  [7] : 

-  initialization: 


di  =  gi  =  e(k) 

W)  (k,  /,  t)  =  ol\  di ,  with  ol\  = 


dfe(k) 


df  R[/j  t\  di. 

-  refinement  iterations  (i  =  2, 3, . . . ,  r)  : 

gi  =  gi-i  +  R [/,  t\  d^_i  Qii-%.,  gradient  vectors 


di  =  gi  +  d^_i  ||gi||2/||gi-i||2,  conjugate  directions 

Oii  =  df  e(k) / (df  R [/,  t\  d0,  step-size. 

We  further  point  out  that  the  rank-r  V-CG  k-f  spectrum  in  (5) 
can  be  interpreted  as  a  Capon’s  k-f  spectrum  (without  matrix 
inversion)  in  a  reduced  r-dimensional  Krylov  space.  That  is, 

p(rank—r )  _  _ f _  /y\ 

eH(k)  D,  (Df  R [/,  t] Dr)-!  Df  e(k) 


ELi|dfe(k)|V(dfR[/,t]d,) 

with  Dr  =  [di,  d2,  •  •  •  ,  dr].  Hence  the  reduced-dimension 
CSDM  becomes  Rd  [/,  t\  =  Df  R[/,  t\  Dr .  Here  the  dimen¬ 
sionality  reduction  matrix  helps  to  alleviate  the  large  eigen- 
spreads  in  the  original  full  dimensional  CSDM  by  re-shaping 
and  subsequently  truncating  the  eigen-values  of  the  original 
CSDM  R [/,  t\,  without  incurring  any  loss  in  signal  detection 
performance  as  long  as  <  Dr  >  captures  the  spatial  modes 
relevant  to  the  steering  vector  e(k).  In  many  applications,  the 
dimensionality  of  subspace  <  Dr  >  can  be  much  less  than  the 
array  dimensionality  N,  and  as  a  result,  adaptive  k-f  analysis 
and  adaptive  beamformers  may  have  greatly  improved  resolu¬ 
tion  performance  due  to  the  re-shaping  and  truncation  on  the 
widely  spread  eigen- values  of  the  estimated  CSDM  R [/,£]. 
Formula  (8)  also  reveals  an  interesting  fact  that  the  rank -r  V- 
CG  k-f  spectrum  can  be  viewed  as  a  weighted  harmonic  aver¬ 
age  of  r  periodogram  k-f  spectra,  each  obtained  by  a  “steering 
vector”  defined  by  a  conjugate  direction  vector. 


2.2.  Data  results 

In  fig.  3(c)  and  fig.4(a),  we  demonstrate  the  high  resolution 
results  from  our  wavenumber- frequency  analysis  using  the  V- 
CG  method((0,  0)  =  (90°,  0°))  with  rank  4  and  5.  Fig.  4(b- 
c)  show  the  k-f  spectrum  plotted  for  (0,0)  =  (90°,  90°)  and 
(0,  0)  =  (0°,  90°).  As  expected,  the  wavenumber  resolution 
is  much  poorer  now  since  we  are  scanning  along  the  cross 
section  dimensions  of  the  array  which  are  small  relative  to 
the  signal  wavelength.  Results  in  figs.  3  and  4  demonstrate 
that  the  conventional  method  suffers  severe  side-lobe  leakage 
from  the  tow-ship  interference  (parallel  modes  in  k-f  slice), 
while  the  Capon’s  method  experiences  mode  smearing  effect 
at  low  frequency  region  and  loss  of  focus  on  spatial  modes 
at  high  frequency  region.  Focused  spatial  modes  can  be  seen 
across  the  frequency  band  of  interest  in  the  proposed  V-CG 
method.  It  appears  that  most  of  the  noise  power  is  primar¬ 
ily  concentrated  in  the  acoustic  region  (marked  by  the  white 
dashed  lines)  with  little  noise  present  in  the  non-acoustic  re¬ 
gion.  This  suggests  that  little  flow  noise  is  present  in  the  array. 

We  also  notice  that  the  slowness  spectrum  provides  a  sim¬ 
ple  way  of  coherently  consolidating  the  power  distribution 
at  different  frequencies  from  the  k-f  spectrum.  For  a  give 
4D  wavenumber- frequency  function  P( k,  /,  t ),  the  slowness 
spectrum  can  be  extracted  as, 

P(W,  f,  t)  =  P( k,  /,  t,  such  that  t?  =  k//).  (9) 

Fig.  5  shows  the  slowness  spectra  of  a  rank-5  V-CG  method. 
Consistent  peak  locations  in  the  slowness  spectra  acoustic  re¬ 
gion  indicate  that  all  acoustic  modes  are  non-dispersive  as  ex¬ 
pected. 

3.  DISCUSSION  AND  CONCLUSION 

Further  work  is  needed  in  effectively  visualizing  the  4-D  k-f 
volume.  Although  “slicing”  the  spectrum  along  each  of  the  x, 
y ,  and  z  wavenumber  axis  is  highly  revealing,  it  is  plausible 
that  some  flow  noise  effects  may  have  been  missed  as  a  result 
of  this  scanning  scheme.  Another  approach  is  based  on  the 
observation  that  because  of  the  cylindrical  symmetry  of  the 
array  (see  fig.  2),  there  should  be  no  “preferred”  orientation 
of  the  wavenumber  spectrum  in  the  y  —  z  plane.  Hence  one 
could  average  out  ky  and  kz  by  integrating  around  a  circle  of 
some  fixed  radius  in  the  y  —  z  plane. 

Finally,  Krylov-based  rank  reduction  V-CG  methods  pro¬ 
vide  significant  improvements  in  performance  over  periodogram 
and  Capon’s  methods  when  applied  to  real  volumetric  data. 
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(a)  Conventional  method  PPEr( k,  /,  t) 


I 


(b)  Capon’s  method  Pcafon( k,  /,  t) 


(c)  Rank-4  VCG  method  (k,  /,  t) 


Fig.  3.  2D  k-f  slices  out  of  4D  wavenumber-frequency  spectrum 
estimates  from  different  methods,  at  a  fixed  time  slot  and  a  fixed 
spatial  direction  ( 0 ,  0)  =  (90°,  0°). 


(a)  PycG*  ^ (k,  /,  t),  viewing  in  kx  direction. 


rank-5  VCG  ky  versus  f 


(b)  PycGk  ^  (k,  /,  t ),  viewing  in  ky  direction. 


(c)  VCG  Pycok  ^  (k,  /,  t),  viewing  in  kz  direction. 

Fig.  4.  2D  k-f  slices  out  of  4D  wavenumber- frequency  spectrum  es¬ 
timate  from  a  rank-5  V-CG  method,  at  fixed  spatial  directions,  (a)  kx 
direction  (0,0)  =  (90°,  0°);  (b)  ky  direction  (0,0)  a  (90°,  90°); 
(c)  kz  direction  (0,  0)  =  (0°,  90°). 


Fig.  5.  Overlaid  slowness  spectra  extracted  from  a  rank-5  VCG  4D 
k-f  spectrum  on  the  3D  volumetric  data,  over  the  frequency  band  of 
interest.  Acoustic  region  is  marked  by  red  dashed  lines. 
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